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We describe "small bodies" in a non-metric gravity theory previously studied by this author. The 
main dynamical field of the theory is a certain triple of two-forms rather than the metric, with only 
the spacetime conformal structure, not metric, being canonically defined. The theory is obtained 
from general relativity (GR) in Plebanski formulation by adding to the action a certain potential. 
Importantly, the modification does not change the number of propagating degrees of freedom as 
compared to GR. We find that "small bodies" move along geodesies of a certain metric that is 
constructed with the help of a new potential function that appears in the matter sector. We then 
use the "small body" results to formulate a prescription for coupling the theory to general stress- 
qq | energy tensor. In its final formulation the theory takes an entirely standard form, with matter 

, propagating in a metric background and only the matter-gravity coupling and the gravitational 

dynamics being modified. This completes the construction of the theory and opens way to an 
f\j ' analysis of its physical predictions. 
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C\| ' I. INTRODUCTION 

<N . 

In general relativity (GR) a test particle moves along a spacetime geodesic. The fact does not need to be added 
as a separate postulate of the theory. Indeed, general relativity possesses a well-defined initial value formulation, so 
to determine the evolution of a body one should just prescribe the initial data for its gravitational field and read off 
{— i ' the trajectory from the resulting spacetime metric. For a "small body" this procedure gives the geodesic motion, and 
the Bianchi identities satisfied by the Einstein tensor are at the root of the derivation. A systematic procedure that 
allows to derive not only the geodesic motion, but also the corrections to it (the so-called gravitational self-force) has 
£SJ , been recently described in 

£> ' In the present author has described a large class of gravity theories that are based on two-forms rather than 
the metric. This class contains general relativity (in Plebanski formulation 0) and can be arrived at rather simply, 
see by dropping the simplicity constraints of Plebanski's theory. A theory from the class is specified by a certain 
"potential" - a scalar function of certain components of the two-form field. Exactly as GR, the theory describes 
just two propagating degrees of freedom. The same class of gravity theories has appeared much earlier in works of 
Bengtsson and Peldan under the name of neighbours of GR, see e.g. ml. These authors' starting point (the pure 
connection formulation [g]) was, however, entirely different from that in J2J, so the equivalence of models proposed in 
[f| to the theory described in (2f is not obvious and was pointed out in [J. 

In the version proposed by this author, the basic dynamical field of the theory is an su(2) Lie algebra-valued two- 
form (complexified in the Lorentzian signature case). With the Lie algebra being three-dimensional, the two- form 
field can be viewed as a triple of two-forms, and these can be declared to span the space of two-forms self-dual with 
respect to some metric. The knowledge of which two-forms are self-dual can be shown Q to determine the conformal 
structure of the metric uniquely. Thus, any theory based on su(2) Lie algebra-valued two-forms is naturally a theory 
of the conformal structure of spacetime. However, it is by no means obvious which metric in the conformal class so 
defined plays a physically distinguished role. Note that in Plebanski formulation of GR this problem does not arise as 
additional simplicity constraints that are imposed on the two-form field guarantee the existence of a preferred metric. 
It is not even clear that there is any physically distinguished metric in the theory. Indeed, this would be a metric in 
which test particles move along geodesies. To find whether there is such a metric, one would need to describe how the 
"usual" matter couples to the gravity theory in question. However, with the theory being that of two-forms rather 
than the metric, this is an unsolved problem. The only case that is currently understood is that of Yang- Mills fields 
which, being classically conformally-invariant, do not help. 

The goal of this paper is to develop the physical interpretation of the theory Q by studying the motion of a 
"small body" . This allows us to sidestep the unsolved problem of coupling to generic matter and develop the physical 
interpretation of the theory remaining entirely within its domain. We shall use the systematic procedure of [If that 
needs only very little adaptation to the theory in question. 

Our main result is that there is a physically distinguished metric in the theory Q along whose geodesies test 
particles move. However, we find that it is the matter itself that supplies the conformal factor that determines this 
metric. Thus, we shall see that the coupling of matter to the gravity theory in question is characterized by a certain 
"mass" function, and the metric in which particles move along geodesies is obtained by choosing the conformal factor 
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so that this function is a constant. If the theory is to preserve the weak equivalence principle the mass function of 
all material bodies must be the same. This requirement introduces a certain universal potential function, see the 
main text. A very similar function appears on the gravitational side, and the theory is thus completely specified by 
prescribing the gravity and matter side potentials. 

Having obtained an expression for the stress-energy "tensor" of a small body, it is not hard to extend it to a 
description of how the general stress-energy tensor of matter couples to our gravity theory. We give such a description, 
thus completing the construction of the theory and making the study of its physical predictions possible. We would 
like to emphasize at the outset that, in spite of the metric appearing in this theory only indirectly, the final formulation 
of the theory is entirely standard: one has usual matter fields moving in a metric background. Only the dynamics of 
gravity, as well as its coupling to matter are modified. However, unlike all previous modification schemes considered 
in the literature, the theory in question modifies both vacuum and non-vacuum GR without adding to it any new 
propagating degrees of freedom. We would like to stress that this feature of the theories considered here is quite 
striking, for it is a rather common belief that the only way to modify Einstein's theory is to add to it new propagating 
modes. 

The organization of this paper is as follows. In the next Section we review the Plebanski formalism for general 
relativity. We also describe how matter (e.g. ideal fluid) can be coupled to gravity in this formulation. In Section 
IIIII we describe the modified gravity theories @ ■ We extend this description to the non- vacuum case in Section IIVI 
In Section [V] we obtain an expression for the stress-energy-momentum two-form of a "small body" and use Bianchi 
identities to determine its motion in Section IVIl An interpretation of the equations we obtain is contained in Section 
IVIII We describe how a general stress-energy tensor is coupled to the gravity theory in question in section IVIIII 
Section ITXl gives a metric formulation of the theory that is most useful for practical applications. We conclude with 
a discussion. 

II. PLEBANSKI FORMALISM 

The aim of this Section is to review the Plebanski formulation of general relativity. As we have already mentioned, 
in this formulation Einstein's gravity becomes a theory of two-forms rather than the metric. Plebanski's formalism 
uses in a deep way the notions of self-duality on two forms (as determined by the spacetime metric) and it is the 
process of "abstracting" this notion from the underlying metric that allows for a deep reformulation of Einstein's 
theory. The original Plebanski's paper [3[ used spinor notations and is not very transparent for a reader who is not 
familiar with spinor techniques. An excellent exposition of the theory is also available in Q, where the problem of 
coupling to matter is discussed as well. 

In this paper, to make it more easy to follow, we will try to avoid using spinors as hard as possible, only resorting 
to spinor techniques when they simplify computations. All such spinor calculations are banned to the Appendix. 

We have also decided to make the exposition of the Plebanski formalism as concrete as possible, so in this section 
we present it as a concrete recipe for deriving Einstein equations. However, before we give such an explicit description, 
let us state the main ideas abstractly. 

Plebanski theory introduces a (complexified in Lorentzian signature) SO (3) vector bundle V over the spacetime M, 
which can be referred to as the self-dual bundle, and a two-form field B l ,i = 1,2,3 taking values in V. The triple 
of two forms B l ,i = 1,2,3 encodes information about the metric g on M via the requirement that B l ,i — 1,2,3 
are self-dual two-forms with respect to g. Indeed, the triple B l spans a 3-dimensional subspace in the space of all 
two-forms, and declaring this to be the subspace of self-dual two forms defines the notion of Hodge duality on two 
forms, which, in turn, can be shown Q to uniquely determine the conformal class of the metric. However, a general 
triple B l of two-forms contains too many components as compared to a metric. Indeed, it needs 3x6 numbers to be 
specified, while a metric has only 10 components. To remedy this, Plebanski imposes the following "metricity" (or 
simplicity) conditions: 

B'ABJ-S^, (1) 

which give 5 equations on the two-form field (the trace of this equation gives the proportionality coefficient and is an 
identity). This brings the number of components in B l down to 13, which is the required 10 components describing 
the metric, plus 3 gauge components related to availability of SO (3) gauge transformations. The volume form of this 
metric is then given by (i/3)(B l A B 1 ). With the two-form field B l being complex, one further needs 10 conditions 
that guarantee that the metric obtained is real Lorentzian. These conditions are given below in ([7]). 

Thus, supplemented by the metricity conditions ([I]) the two-form field contains just the right amount of information 
to describe a metric. One now needs a second order differential equation on B % . To obtain this one notices that there 
is a unique connection A 1 satisfying: 



D A B l = 0, 



(2) 
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where D A B' := dB l + e ijk A j A B k . Indeed, this gives 4x3 algebraic equations for 4x3 components of A 1 , which 
fixes it uniquely, provided a certain non-degeneracy conditions for B l are satisfied. We shall refer to this connection 
as incompatible and denote it by As- When B % satisfies the metricity conditions fl} the ^-compatible connection 
turns out to be equal to the self-dual part of the metric-compatible one. One can now compute the curvature 
F l := dA % + (l/2)e l: > k A> A A k of the B-compatible connection. A natural second order equation on B l is obtained 
by requiring that the curvature F l is "proportional" to the two- form field B l . However, it is natural to allow for the 
"proportionality coefficient" to be an "internal" tensor: 

/•" <!>'■' /,". (3) 

where the quantities <5 y are at this stage arbitrary. It can then be shown that the "internal" tensor <E> y must be 
symmetric and that its trace part must be a constant: 

= + hs ij A, 9 ij = <F fo) , =0, A= const. (4) 

o 

The relation ([3]) then becomes a set of 18 equations for 13 components of the i?-field as well as for 5 (traceless 
symmetric) undetermined components of . To see that ^ is equivalent to Einstein equations one notices that 
it states that the curvature of the self-dual part of the spin connection is self-dual as a two-form. This is known to be 
equivalent to the Einstein condition. 

Finally, let us note that all equations of Plebanski theory can be obtained as Euler-Lagrange equations for the 
following action: 



S[B, A,*] = J B l A F\A) - i ^ + A B l A B j . 



(5) 



Indeed, the variation with respect to the traceless tensor ty 13 ' gives (JTJ , variation with respect to the connection gives 
(|2"|) , while variation with respect to the two- form field gives the main dynamical equation ([3]) . 

To make the above description little less abstract let us reformulate it as a concrete recipe for writing down Einstein 
equations once a spacetime metric is given. The starting point of Plebanski method of deriving Einstein equations 
is the same as in the tetrad method: one has to find a suitable tetrad. For a diagonal metric there is no ambiguity, 
but for a non-diagonal one it is possible to use the available freedom of Lorentz rotations to bring the tetrad to a 
convenient form. Thus, we assume that we have found a convenient collection of one- forms e 7 ,i = 0, 1,2,3 so that 
the metric is: 

ds 2 = e 1 ®e J r)u, (6) 

where r\u is the Minkowski metric. 

The second step is to form a set of three two-forms B l ,i — 1,2,3 which are self-dual with respect to the given 
metric, and satisfy: 

B l AB^S ij , B i A (B° )* = 0, Re(B i A B i ) = 0, (7) 

where (B 1 )* are the complex conjugate two-forms. This task is easy if one has a tetrad in one's disposal, with a 
possible solution being: 

B 1 = ie° Ae 1 ^ 2 A e 3 , B 2 — ie° A e 2 — e 3 A e 1 , B 3 = ie° A e 3 — e 1 A e 2 . (8) 

It is easy to see that all the required conditions ([7]) are satisfied, and that the above two-forms are indeed self-dual 
with respect to ([1}. It can also be shown that the converse is true: given a triple of two-forms B % satisfying there 
is a unique real metric with respect to which the two-forms B l are self-dual. 

The third step is to find an su(2) connection A 1 that is "compatible" with the above set B % of two-forms, in the 
sense that the covariant derivative of B % with respect to A 1 is zero: DaB 1 = dB l + e 1 ^ A* A B k = 0. To obtain such a 
connection one has to solve the following system of linear algebraic equations for the components of the connection: 

dB 1 + A 2 A B 3 - A 3 A B 2 = 0, dB 2 + A 3 A B 1 - A 1 A B 3 = 0, dB 3 + A 1 A B 2 - A 2 A B 1 = 0. (9) 

In doing so we first have to take the exterior derivative of the three two forms B l , and then write down the above 
system of 3 x 4 equations for the connection components a\ : A 1 — a^O 1 . This is an exercise in algebra, of roughly 
the same degree of complexity as one that arises in the determination of the rotation coefficients in the tetrad-based 
approach. What simplifies the game somewhat is that in the tetrad-based approach, at least in principle, one has 
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4x6 equations to write, for the same number of the rotation coefficients, while in the case of Plebanski formulation 
there is just half that number. No information is lost, however, as all quantities in the Plebanski case are complex. 
In practice, for a given metric (usually possessing some symmetry properties), most of the connection coefficients are 
zero by symmetries both in the tetrad and the two- form cases, so the amount of work one has to do to find A 1 is only 
a little less than in the tetrad-based scheme. As in the tetrad-based scheme it is much easier to verify a conjectural 
solution than find one, as the latter involves some guesswork on which of the components are zero. Finding the 
connection is the most laborious part of the computation. 

The fourth step is to compute the curvature two-form F l = dA % + (1/2) e y A-? A A k . In components: 

F 1 = dA 1 + A 2 A A 3 , F 2 = dA 2 +A 3 AA\ F 3 = dA 3 + A 1 A A 2 . (10) 

This is a simple exercise in differentiation. The step of computing F l should be compared to the curvature computation 
in the tetrad-based approach. In that case one needs to compute six two-forms, as compared to only three in the 
Plebanski case. 

The fifth step is to write all the two-forms that appear in F l in terms of the six basic two-forms that are used in 
([8]). This may involve some algebra in the case the metric is non-diagonal. This step is exactly the same as in the 
tetrad-based scheme, when one writes the curvature components in terms of the basic tetrad two-forms (before the 
Ricci tensor can be found by contracting a pair of indices). 

The final step is to replace the tetrad two-forms appearing in the result for F l by their expressions in terms of B i 
and (B 1 )* . It turns out to be more convenient to use the two- forms B % — —(B 1 )*. We have: 



eW^l^+B 1 ), 
2i 

e Q Ae 2 = ^(B 2 + B 2 ), 
2i 

e°Ae 3 = -(B 3 + B 2 ). 
2i 



e 2 A e 3 = 


V 

2 V 




e 3 A e 1 = 


-(B 2 


~B 2 l 


e 1 A e 2 = 


-(B 3 
2 y 


— B 3 ). 



(11) 



This last step does not have a direct analog in the tetrad-based method. 

This is it! Once an expression for the curvature F l in terms of the self- and anti-self dual two- forms is found, one 
can immediately write down the Einstein equations. Indeed, we have obtained the curvature in the form: 

F i = M ij B j + N ij B j , (12) 

where M lJ , iV u are some matrices built from the components of the metric and their first and second derivatives. The 
Plebanski field equations read: 

F* = (V J + ^'Aj B J + 4nGT\ (13) 

where is a traceless matrix, A is the cosmological constant and T l is the "stress-energy-momentum two-form" , 
see below, zero in vacuum. Thus, the vacuum Einstein equations are simply: 

TrM = A, N lj = 0, (14) 

which gives ten equations, as it should. As a bonus, we also obtain the components of the (self-dual part of the) Weyl 
curvature tensor. Indeed, we have ^ = (M 1 -?)^), with (tf) standing for the tracefree part. 

Let us also write the general non-vacuum equations. To this end we need a general expression for the "stress- 
energy-momentum two-form" in terms of the quantities characterizing the gravitating matter. A general Lie-algebra 
valued two-form admits an expansion of the type (TT2"|) , thus giving rise to 9 + 9 = 18 components. However, in 
general relativity the "right-hand-side" of Einstein equations - the matter stress-energy-momentum tensor - has ten 
components. Thus, the stress-energy-momentum two form of Plebanski formalism cannot be a general Lie-algebra 
valued two-form. It needs to satisfy: 

/ ' ' />" - <)''. (15) 

which should be compared to the first of the conditions in ([7]), plus certain reality conditions, see below. The conditions 
(TTSJ) imply that T % has the following structure: 

T l = -TB l --T ij B j , (16) 
6 2 K ' 
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where the numerical prefactors are added for future convenience, and T, T i] are an arbitrary scalar and 3x3 "internal" 
tensor correspondingly. This gives in total 10 components in T l as should be expected from a quantity representing 
the usual stress-energy tensor T M „. The interpretation of the components in (fT6| is as follows. The scalar T is just 
the trace T = g^T^ of the usual stress-energy tensor. The 3x3 tensor T lJ can in turn be decomposed into its 
symmetric and anti-symmetric parts: 

T ij =a ij +ie ijk u k , (17) 

and takes care of the traceless T M „ — (l/Ajg^T part of T M „. The symmetric part a 13 turns out to have the meaning 
of the spatial components of the stress-energy tensor, i.e. characterize the stress of matter. The anti-symmetric part 
u % gives the momentum vector. The quantities T,a 13 ,u l are all required to be real. Thus, the symmetric part of T %3 
is required to be real while its anti-symmetric part is purely imaginary. 
The non-vacuum Einstein equations then take the form: 

TrM = A + 2itGT, N 13 = -2%GT ij . (18) 

For many applications one is interested in the simplest type of matter - that given by the perfect fluid. For the 
perfect fluid only the trace part of the stress tensor is different from zero: a' 3 ~ S' 3 . It is thus completely characterized 
by its energy p and pressure P densities and the momentum vector it\ The corresponding stress-energy-momentum 
two- form is: 

T} luid = \{P~ 3P)B l ~\{{P+ P)S ij + ie ijk u k ) &. (19) 



III. MODIFIED GRAVITY: THE VACUUM CASE 

We will start our presentation of the theory 0] by describing the vacuum case. We follow closely a recent description 

Dt- 

The class of theories in question can be obtained by relaxing the simplicity conditions (fl]) of the Plebanski theory. 
Thus, the main idea is to allow all the components of the two-form field B % to become dynamical. Recall that a general 
two- form field B % determines a conformal structure of spacetime by requiring that the triple B l , i = 1, 2, 3 is self-dual 
in this conformal structure. Thus, we are going to obtain a theory in which the spacetime conformal structure, not 
the metric becomes the main dynamical object. 

How can one get natural field equations describing the dynamics of B l l As in Plebanski formulation of GR one first 
computes the B-compatible connection A B and then its curvature F(A B ). As before, it is natural and instructive to 
decompose the curvature into the basis of two- forms B l and B % ; 

F\A B ) = M ij B j + N i3 B 3 . (20) 

The most natural field equations are the same as in the Plebanski case Q. Thus, we require the curvature of the 
connection As to be purely self-dual: 

F l {A B ) = <f> ij B j <^> M ij = & j , N ij = 0, (21) 

where <fr lJ is some purely gravitational tensor to be described below. The system of equations (|21[) gives us 18 equations 
for the 18 components of the two- form field B % . However, it also contains the so-far unspecified functions $ y and so 
is not complete. Note that in the GR case we have exactly the same system of 18 equations, but in that case for 18-5 
quantities B l (the two-form field B % modulo the conditions B l A B 3 ~ S 13 ) . In addition the trace part of $ u is either 
zero (no cosmological constant case) or constant (the cosmological constant), and is thus not an unknown field. The 
system of 18 equations is thus that for 13 components of B l and the remaining 5 components of . 

In the general case the system of equations (|2Tj) can be completed by considering analogs of "Bianchi" identities. 
Thus, we note that the components of AP 3 , N %3 in (|2U)l are not independent. Indeed, we have the following Bianchi 
identity: 

D Ab F{A b ) = 0. (22) 

This gives: 



{D Ab M 13 ) A B 3 + {D Ab N* 3 B 3 ) = 0. 



(23) 
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Another important identity is obtained by using the compatibility equations Da b B 1 — 0. Taking another covariant 
derivative and using the definition of the curvature we get: 

e iJ k F j (A B ) A B k = <^=> e l ^ k M jl B l A B k = 0. (24) 

This last equation can be conveniently interpreted as follows. Let us define a conformal "internal" metric: 

I!' ■ I f //'•''. (25) 

Then (|24|) can be rewritten as: 

e ijk M jl h lk _ q (26) 



Let us now also introduce an action principle that leads to (|21[) as Euler-Lagrange equations. This is easy to write, 
we have: 

S[B, A, $] = J B i A F'(A) - i^y'B' 1 A (27) 

Varying this with respect to A 1 we get DaB 1 — 0, which allows to solve for A in terms of B, varying the action 
with respect to B % we get (|21[) . We also note that only the symmetric part of the field 3> y enters the action, so it is 
necessary to assume that <& y in (f2Tj) is symmetric. 

It remains to clarify the meaning of the variation with respect to To these end we shall use the Bianchi 

identities (f2"3")) . (f2"4")) . Using (j2"3")l and field equations (f2"Tj) we see that we must have: 

Di B $' 3 AB 3 ' = 0. (28) 

Let us multiply this equation by the one-form l^B 1 and sum over i. Here £ is an arbitrary vector field and l^B 1 is 
one-form with components (i^B 1 )^ := ^ a B l afl . However, for any vector field £ we have: 

A = h^B 1 A B j ) ~ h", (29) 

where h lJ is the internal metric introduced above. This gives us the following equation: 

h ij D AB $ ij = 0. (30) 

Now, using the symmetry of we can rewrite this equation as /i y + 2e lkl A k ^ ) — 0. However, the other 
Bianchi identity (f2"4")l together with the field equation M y = $ y implies e $ ^hP % — and so we must have: 

Wd&i =0. (31) 

The identity (|3"Tj) implies that the quantities and <j> y are not independent. This can be seen quite clearly by 
considering the last term in the action (|27p. Using the introduced above tensor h lJ we can write the integrand as 
V := h^^i times some volume form. We then get: 

dV = <!•"<//," + //",/<!>' ' = <f> ij dh ij , (32) 

where we have used pip . This means that (i) the last term in the action is only a function of the h 1 ^ components of 
the two-form field B l ; (ii) the quantities <J>* J ' are expressible through and are given by: 

*« = (33) 

Below we shall characterize the "potential" V(/i) in more details. For now let us note that having expressed the 
unknown functions $ y in terms of the the components of the two-form field B % we have closed the system of equations 
(l2"Tj) . as it is now a system of 18 equations for 18 unknowns - components of the B % field. 

To understand the structure of the potential V(h) it is convenient to parametrize the "internal" metric h lJ by its 
trace and the traceless part: 

h ij = ^Tr(h) (S ij +H ij ) , (34) 
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where if 1 - 7 is tracefree. It is then easy to see that for any function = /(Tr(/i), if 1 - 7 ') 

9/ / 9/ 9/ if fci \ , 3 9/ 



d/i^ \9Tr(/i) dH kl Tn(h)J Tr(h)dH 

In particular, we have 



** + ^ a ^,- (35) 



9/ hij = t^ttt 1 ^)- ( 36 ) 



dTi(h) 
Thus, one has: 

V = ^=Tr { k)^ y (37) 

where is given by (|33p . Thus, we learn that the potential must be a homogeneous function of order one in its 
argument Tr(h): 

V(h) = A^-U(H), (38) 

where U is a dimensionless function that only depends on the tracefree part H lJ of h %3 , and is normalized so that 
U(0) = 0, that is: 

U(H) = l + |Tr(if 2 ) + 0(if 3 ), (39) 

where a is some dimensionless parameter. The quantity A is a constant of dimensions 1/L 2 that needs to be introduced 
to give V(h) the correct dimensions. Below it will be identified with the cosmological constant. Using our definition 
of the "internal" metric h v we can write: 

B* = 7^(B k AB k ), (40) 



which of course defines /V-? only up to a conformal factor. One can now write down the action (|27|) as a functional of 
only the two-form and the connection fields: 

S[B, A] = J B i A F\A) - hW^-B 1 A B\ (41) 

where if 4 - 7 is defined as the traceless part of the internal metric /V- 7 , and is independent of the conformal freedom 
present in the definition of h % K Note that this action is an off-shell one, that is it can be varied with respect to the 
dynamical fields B l , A 1 to obtain field equations. Alternatively, one can work with a version that uses extra "Lagrange 
multiplier" fields, see below. 

Note that the dimensionfull constant A here should be identified with the cosmological constant of our theory. 
Indeed, one can, e.g., consider the metric describing a homogeneous isotropic Universe. In such a Universe if = by 
symmetries, and so the metric evolves exactly like in general relativity with the cosmological constant A. In other 
words, a solution of our theory describing a homogeneous isotropic Universe is the same as in GR with cosmological 
constant A. Thus, to define the modified gravity theory in question one only needs to specify a dimensionless function 
U (if) of a dimensionless traceless tensor H lJ . All physical dimensionfull parameters present in the theory are as in 
general relativity. Let us remark that we could have included A into the definition of the potential U(H) and thus 
made it dimensionfull. However, as we shall see below from matter coupling considerations, it is more convenient to 
make the potential function dimensionless, for one can then use similar potential functions in both the gravity and 
the matter sectors. 

In terms of the introduced potential U(H) the main set (|2ip of field equations becomes: 

F i {A B ) = A f -^^7 + ~S ij lf\ B j , (42) 



dH^ 3 

where we have introduced the Legendre transform U of the potential U 



(43) 
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The function U can be viewed as either that of or of the quantities: 

dU 

**J/A := ^rr. (44) 

The theory with an arbitrary "cosmological function" A(^) := AC/(W/A) defined by the Lagrangian 

S[B, A] = J B l A F i {A) - ~ - ^Af*)^ B i A B\ (45) 

is that of the original paper 0] • Field equations (|42p are most easily obtained precisely in this "Lagrange multiplier" 
formulation. However, the viewpoint suggested by (|4"Tj) . namely that of the gravity theory being the BF theory (the 
first term in (|41|) ) plus a potential term for the H lJ components of the two-form field will be more convenient for our 
purposes here. 

To summarize, we have seen that the condition B % A ~ 6^ of Plebanski formulation of GR can be relaxed and 
how the Bianchi identities still lead (in a unique way) to a consistent theory. Note that what one obtains is a class 
of gravity theories rather than one theory, for a theory is now specified by a choice of the dimensionless "potential" 
function U(H) of the components H v of the two-form field B l . The potential can be completely arbitrary. One can 
obtain back general relativity (in Plebanski formulation) by making the potential function U(H) infinitely steep so 
that the quantities iP J are set to zero. However, if one sets U (H) = 1 one obtains a topological theory - the so-called 
BF theory with a cosmological constant. 

Having achieved a formulation of the theory in vacuum it is very important to continue to develop the theory and 
allow for a non-trivial right hand side of our equations - for matter to be present. Indeed, pure gravity is only of 
academic interest, and the real world gravity is both produced and felt by material bodies. 



IV. MODIFIED GRAVITY: NON- VACUUM CASE 



The vacuum field equations (pH]) of modified gravity were exactly the same as those of the vacuum Plebanski 
theory. As in Plebanski case, it is natural to describe the effect of matter on the modified gravity "geometry" by 
allowing a non-zero T l to be present on the right hand side of (|2ip . This is in the spirit of Einstein equations, where 
the stress-energy-momentum of matter appears on the right hand side of an equation for the curvature and thus affects 
the geometry. This in turn implies that the stress-energy-momentum of matter should satisfy some conservation laws 
which to a large effect determine its motion in a given background. 

Thus, we shall keep the field equations ([3]) as our main dynamical equations even in the case of non-zero T % . 
However, now that we have removed the condition B % A B? ~ it no longer consistent to impose the condition 
T l A B J ~ 8 V] either. In general, the matter "stress-energy-momentum" two-form T l will have all components: 

T* = ^Q ij B j T ij B j , (46) 

where the interpretation of components T r] is similar to that in GR, see (jTTJ) , and Q lJ are some new "internal" 
components of T\ We note that, in the case Q y ~ <5 y , the trace T of the usual GR stress-energy-momentum tensor 
is just a multiple of the trace Tr(Q). The quantity Q tJ , by analogy with T u can be referred to as the "internal" 
stress-momentum of matter. Let us decompose Q lJ into its symmetric and anti-symmetric parts: 

QV = x i i +ie ijk £ k . (47) 

The quantities x tJ ? C then receive the interpretation of "internal" stress and momentum correspondingly. However, 
at this stage, there are no reasons to require x %3 ; C to be real, while the similar quantities in the "spacetime" stress- 
momentum T u are observable quantities and thus must be real. We note that the various pieces in the decomposition 
(I46p have the interpretations of stress and momentum (internal and "spacetime" ones), but there is no energy density 
anymore, the later appearing as the combination of the traces of the internal and spacetime stress tensors. For this 
reason, and also for brevity, we shall refer to T % as the stress-momentum two-form from now on. 

Let us now rewrite our field equations as relations between the curvature and stress-momentum components: 

F\Ab) = $ ij B j + AirGT* M ij = & j + 2t:GQ 1 \ N lJ = -2irGT l \ (48) 

where as before the "gravitational" quantity ^ is at this stage arbitrary and is to be determined via the help of 
Bianchi identities. In the field equations (|4"5)) the quantity G is the usual Newton's constant. Note then that the 
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quantities Q u , T u must have the dimensions of energy density, so the above interpretation of the components of Q y 
as the "internal" stress- momentum is consistent with their dimensions. 
Now, let us, as before, construct an action that leads to (l48l). We have: 



S[B, A, $, ro ] = J B i A F\A) - A B ] - 4wGS m [B, m ], (49) 

where (j) m is a collective notation for all the matter fields, and S m [B, <j) m ] stands for the matter part of the action, 
which is assumed to depend on the "gravitational" background only via the two- form field B % . In principle, one 
can also envisage the possibility of the matter fields (e.g. fermions) coupling directly to the connection A 1 , but this 
does not bring anything conceptually new, only complicates the analysis, so we shall not consider this possibility any 
further. Defining: 

- f| (50) 

we get the equation (|48[) when varying the action with respect to the two-form field B % , and the compatibility equations 
D A B l — when varying the action with respect to the connection. 

The next step in interpreting the above theory is to note that the stress-energy-momentum two-form T l satisfies 
some conservation laws. Indeed, since the matter part S m of the action must be diffeomorphism invariant, the following 
identity must hold: 

= 6 6 S m = J ^h'Pm + J T l A 8 t B\ (51) 

where 8^ is a variation of fields under an infinitesimal diffeomorphism generated by a vector field £. The first term in 
(|5ip vanishes by matter equations of motion, while the other term gives: 

= jT l AD AB LtB\ (52) 

where (,£ denotes contraction of a form with a vector field, so e.g. i^B 1 is a one-form with components (l^B 1 )^ := 
£ a B l a ^. This expression follows from the following formula for the action of diffeomorphisms on B % : 

C^B 1 = l ( D a B 1 + D a l ( B\ (53) 

where A can be taken to be arbitrary, and the fact that D A B l = for A — As- Now integrating (|52[) by parts and 
taking into account that £ may be of compact support, we can conclude that the integrand must vanish: 

D Ab T 1 A l ( B 1 = 0. (54) 

This should hold for any vector field £, so we get four "conservation" equations. 

Another important identity that we can obtain for T l follows from gauge invariance of the action. Thus, we similarly 
write (IBTj) but now consider the variations of the fields under a gauge transformation. For the B l field this is: 

8,ir (55) 

where lu ij is an infinitesimal anti-symmetric matrix - a Lie-algebra element of SO(3). We can therefore conclude that 

rAw¥=0 (56) 

for any matrix u> and thus: 

A B^ = 0. (57) 
In the language of the decomposition (146[) this translates into: 

QW h m = o ^ ,'- !i '(J !l h lh = 0. (58) 



This identity is satisfied by any T % that follows via ([50)) from a gauge-invariant action, and will be of importance 
below. 
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We can now derive a Bianchi identity for the gravitational quantities $ y . Let us rewrite (|23j) as: 

D AB ¥ j A B j + AttGD Ab T 1 = 0. (59) 

Let us now take the wedge product of this expression with the one- form i^B 1 . Using (|54l) we see that the field $ w 
must satisfy: 

L(.B l A D Ab $ y A B j ~ h ij D Ab = 0, (60) 

where we have used (|29|) . We can now write down all the terms in the expression for the covariant derivative D Ab 
Similarly to what we had in the pure gravity case, there is a term containing the connection A l B and proportional to 
& l h . However, we can again conclude that this term is zero. Indeed, the second Bianchi identity (j2"4")l together 
with field equations (|48|) says: 

e «fc($ii + 2nGQ jl )h lk = 0. (61) 



However, we have seen above that the invariance of the material action under gauge transformations implies J58J), and 
thus the second term here is zero, which implies that the first term is zero as well. Thus, from (|60p we conclude that: 

h ij d<f> ij = 0, (62) 

which is exactly what we had in the pure gravity case. All the remaining steps from the previous Section go unchanged: 
we arrive at conclusion that the potential term in the action proportional to V = h 10 ^ is a homogeneous function 
of order one in Tr(/i) times the cosmological constant A times an arbitrary dimensionless function U(H). Thus, the 
full theory is obtained as simply the gravitational plus matter parts: 

S[B,A,<f> m ] = J B i hF i {A)-KWlB i hB i -^GS m [B,(t> m ]. (63) 

This solves the problem of coupling of the class of generalized gravity theories we have been considering to generalized 
matter. It only remains to supplement the matter part of the action with some appropriate reality conditions, for it 
is in general complex. A way to do this is to require the components of T l that are "directly" observable to be real. 
However, to be able to do physics with our gravity theory we need to understand how material actions S m [B, <j) m ] can 
be formed and which stress-momentum two-forms T l can arise in our non-metric theory. 



V. STRESS-MOMENTUM TWO-FORM OF A "SMALL" BODY 

In the previous Section we have seen that the matter stress-momentum two-form T 1 satisfies the following "conser- 
vation" equations: 

L^B i A D Ab T 1 = 0, T [l Afl j] =0. (64) 

In GR the stress-momentum has the special form (fT6"| and the second of these equations is automatically satisfied, 
while the first gives the usual conservation of energy equation when the two-form field B l is metric B l A B^ ~ <5 y . 
The conservation equation can then be used to conclude that "small bodies" move in GR along geodesies. A very 
important question for the theory just developed is what the notion of geodesic generalizes to in the case of an 
arbitrary two-form field background. To understand this we shall employ the same methods as are used in GR. 

A particularly efficient method that allows to study this question has appeared recently in a paper by Gralla and 
Wald [l[ . This work employs the machinery of asymptotic expansions to derive results on motion of "small" bodies 
in GR both in the leading approximation, which gives the result that bodies move along geodesies, as well as in the 
sub-leading one, which leads to results on gravitational self-force. For our purposes we only need the analog of the 
first of these. Thus, most of the machinery developed in this work is actually unnecessary here. However, some key 
ideas of Section IV of this paper will still be used. 

As in [l[, the first step is to derive, using the asymptotic expansion techniques, that the stress-energy- momentum 
tensor has a well-defined limit approaching a distribution. To this end, consider a family B l (\, x a ) of two-form field 
backgrounds. Here A > parametrizes members of the family and x a are coordinates of some convenient coordinate 
system. We shall assume there exist coordinates such that B l (X, x a ) are smooth in both A and x a at least sufficiently 
far away r > RX,r 2 — J2i=i( xl ) 2 from the particle, where R is some universal constant, and that for all A, r > R\ 
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the two-form B l (X, x a ) is a solution of vacuum field equations of our theory. Let us consider the following expansion 
for the two- form field as A — * 0: 

B i {X) = B i + Xb i + 0{\ 2 ), (65) 

where for now both the background B % and the "perturbation" b % are functions of all coordinates x a , the background 
B l satisfies the vacuum field equations F 1 (Ab) = ^ B 3 , and the "perturbation" satisfies the linearized field equation, 
at least sufficiently far from the particle r > e. 

Now given a background B % satisfying the field equations, let us define an operator G*(&) via: 

F\A B+Xb ) - <!>''://' + XV) = XG\b) + 0(X 2 ). (66) 

In terms of G 1 the linearized field equations read G l (b) = 0. 

By definition, the stress-momentum two-form is defined as a distribution on spacetime whose action on an arbitrary 
smooth Lie-algebra valued anti-symmetric tensor f 1 is given by integrating the right-hand-side of linearized field 
equations, i.e. G l (b) against /\ Or, using the fact that G l is sclf-adjoint, we can define: 



4ttGT(/) = / G l (f)Ab\ (67) 
Jm 

where b % is as in (|65|) . 

The definition of the notion of self-adjoint that is used here is as follows. It can be shown that for arbitrary 
two- forms b l ,f l the four-form: 

b i AG i (f)-G i (b)Af=dX(b,f) (68) 

is a total derivative, where X(b, /) is a certain 3-form depending on both b and /, as well as the background B l . We 
shall not attempt demonstrate this property in the present paper, as the computation is quite technical and it would 
take us too far. 

If we now integrate the expression (|68|) over a region r > e and use the fact that b % satisfies the linearized field 
equations G l (b) = 0, r > e. We get: 



Taking the limit e — > we get: 



b l AG\f)= / X(bJ). (69) 

r>e J r—e 

Similar to what happens in the GR case, it can be seen that this limit exists and is different from zero if b l ~ l/r as 
r — ► 0. Thus, when this is the case, the stress-momentum distribution can be written as: 

T(/) = Jdte^T^(t)r pa (t). (71) 

Here T l (t) is a Lie-algebra valued two- form along the curve 7 (given in the chosen coordinate system by x % = 0), and 
p(t) := f l (t, r = 0) is the value of the test two-form f l along the curve 7. We note that to write this expression we 
have chosen a particular background metric (see below on how this is done). Also note, for future use, that under 
conformal transformations of the background metric dt — * Vtdt,e^ pa — > f2~ 4 e MI/pcr , and so for the above distribution 
to be invariant under conformal transformations of the auxiliary background metric the stress-momentum two-form 
must transform as (t) — > 51 3 T^ l/ (t). We shall use this fact in the next Section when we check the behaviour of the 
evolution equations under conformal transformations. 

Now, using the Bianchi identity that holds in our theory, one concludes that the above distribution must vanish on 
test two-forms f l of the form: 

D Ab l s B\ (72) 

where £ is an arbitrary vector field and B l is the background two- form field appearing in (|65|) . As in [l|, to derive 
consequences of the arising "conservation" equations, we shall first consider the case of special vector fields of the 
form: 

e =x a F(x 1 ,x 2 ,x 3 )c> 1 (t), i = 1,2,3, (73) 
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where F(x 1 ,x 2 ,x 3 ) is an arbitrary function such that F(r = 0) = 1. As we have already mentioned the coordinates 
x l are chosen in such a way that the curve 7 along which the body is moving (i.e. in the neighbourhood of which 
vacuum field equations are assumed) corresponds to r = 0, where as usual r 2 = Ysi=i( xt ) 2 - This still, however, leaves 
a considerable freedom in the choice of the coordinates. Let us use the background two-form field B l to help with 
this. Thus, recall that B l defines a conformal metric. Since any metric is locally flat, as in [l|, we can always choose 
the coordinates locally so that this conformal metric is just the Minkowski metric. This means that, without loss of 
generality we can assume the two-form field B % in the small neighbourhood of 7 to be given by: 

B l = A*Bi (74) 

where the two-forms B- are those describing the Minkowski spacetime: 

= idt A dx z - - ^eH k dxi A da£, (75) 

and are arbitrary matrix- valued functions of spacetime coordinates. Note that we have introduced a new type of 
indices - underlined ones, to distinguish between the "internal" SO (3) bundle where the original fields take values and 
the "metric" bundle where the metric two-form field (]T5[) lives. The covariant derivative Da only acts on the original 
non-underlined indices. The matrix A" is defined modulo conformal rescalings of the metric introduced, which sends 
B l m — ► SI 2 B^, an( ij since the background two- form field is by itself independent of any choice of the metric, transforms 
A" — * f2^ 2 A". One can also do a Lorentz rotation on the metric two-forms ([75]) that acts on the underlined indices. 
Thus, the quantities A" are only defined modulo such conformal rescalings and SO (3) rotations. 

Now using the special vector fields ([73]) in the test two- form (|72|) . and taking into account that only the term in 
which the exterior derivative acts on the coordinate functions x l gives a non-zero contribution in the limit r — > 0, we 
get that the stress-momentum distribution must vanish on the following set of test two-forms: 

A"-(t)dx L A {ic°{t)dx l - - i£{i)dt + eH k d{t)dx k ) (76) 

for any choice of Z and functions c°(t),c L (t). This gives us 3 x 4 conditions on T l (t) that we would like to exploit 
to deduce the form of the stress-momentum. Here and in what follows the notation f(t) stands for the value of the 
function / along the curve r = 0. Thus, f(t) := f(t, r — 0). 

It is now convenient to consider a related quantity: T L (t) := T z (t)A"(t). The equation in question then becomes: 

T l -(t) A dx l - A {ic°(t)dxi - id{t)dt ~ S- k cL{t)dx k ) = 0. (77) 

Let us decompose: 

T l -(t) = iA^-(t)dt A dxi - lB i i(t)ei M -dx k A dx 1 -, (78) 

where A^(t),B^(t) are some unknown matrix- valued functions of time. Setting c-{t) — and thus extracting the c° 
component of the conservation equations we immediately get: 

A& ] = 0. (79) 

Setting c°(t) = and c-(t) ~ 6— we get, after some algebra: 

B— + Tr{A)S 1 -^ - A— = 0, (80) 

where we have suppressed the dependence on t for brevity. From this equation we immediately conclude that B lJ is 
also a symmetric matrix, and that Tr(£?) = — 2Tr(^4), while the traceless parts of A^^B 1 ^- are equal. Let us denote 
these traceless parts by x^(i), and (a multiple of) the trace of say B 1 ^ by m(t). Then we obtain the following form 
of the stress-momentum distribution: 

T%) = \m{t)B± - \m{t)& + (81) 
D 2 

where we have used the definition ([75)1 of the Minkowski space two-forms, and B 1 := — (B 1 )* is the anti-self-dual 
two-forms. It is instructive to compare this result to the GR one. In that case A 1 - = <5 4i , and no non-trivial self-dual 
part of the stress-momentum tensor is possible, so x^- — 0- The remaining two-form is that corresponding to the 
ideal pressureless fluid, see (flUj) . with coordinates chosen such that the momentum u 1 = 0, as it should. We have thus 
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recovered the GR result, formula (45) of [1(. We see that the main modification arising in our case is the presence of 
an arbitrary traceless part x v - m the self-dual part of the stress-momentum two-form. 

The second "conservation" equation in (|64|) can also be exploited. After some simple algebra we find that along 
the curve r = it is equivalent to the condition: 

e yfc(A- 1 )^(t)x^(i)A J 'i(i) = 0. (82) 

Let us also note the transformation properties of the quantities that appeared in (jSTj) . Since T % — > f2 3 T 4 , we have 
m(t) — ► CI m(t), x^(t) — > f2 _1 ;<^(i), which are the correct transformation properties for the quantities having the 
dimensions of mass. 



VI. MOTION OF A "SMALL BODY" 



Having extracted the form (|8ip of the stress-momentum two-form T 1 of a "small body" we are ready to find equations 
that such a body must satisfy during its motion. Thus, we are looking for an analog of the GR statement that "small 
bodies" move along geodesies. To this end we once again use the fact that the stress-momentum distribution, whose 
form (|81[) we have determined above, must vanish on test two-forms of the form Da b l^B 1 , where £ is an arbitrary 
vector field. 

The computation one has to do is conceptually clear, but a bit involved. A particularly efficient way to do it is to 
use spinors. However, the resulting intermediate formulae are not particularly transparent for somebody not familiar 
with spinor techniques. For this reason we shall use a shortcut based on the fact that changes to the final result only 
come from the self-dual part of the stress-momentum two-form, and the contribution of the anti-self-dual sector is 
completely unchanged from the GR case. This has to be verified, and, as we have said, the easiest way to do this is to 
use spinors. We give a complete derivation in the Appendix. Here we only deal with (the most interesting) self-dual 
part, which can be easily done without spinors. 

Thus, we decompose the stress-energy distribution two-form T 1 into its self- and anti-sclf-dual parts and write: 

i dt(Ti d {t)+Ti sd {t)) ADa^B* = 0. (83) 

This must hold for any vector field £. 

Let us analyze the T l sd {t) term of (|53")) . In the previous subsection we have found that this part of T l can be written 
as: 

TUV) = ^Q ij B j , (84) 
where B l is the background two-form field and Q lJ is a tensor given by: 

Q ij = (A-JHt) (~m(t)5V- + y*l(t)\ (A" 1 )", (85) 

where X^-(t) is as introduced above in (|8ip and the matrix (A -1 )- 7 ^ on the far right is not evaluated at r = 0. We will 
only be interested in this tensor along the curve r — 0. Then the quantity Q*- 7 (t) is symmetric Q^(t) = because 
X^it) is symmetric. We also get the following simple expression for the "mass" m(t): 

K l i{t)Q l i{t) = m{t), (86) 

where h lJ is the "internal" metric now defined as: 

h U = A ik A jk. ( 87 ) 

We also see that the second "conservation" equation (|82|) becomes in terms of Q y : 

e ijk Q jl (t)h lk (t) = 0. 



Now using the fact that Da b B % = we can take the two-forms B J in (|84[) under the operator of covariant derivative. 
We then use the identity (|29[) to get for this first term: 



m(t) 



\ J dtQ ij (t)D AB ^(B i A B j ) ~ j 'dt (lQ i i(t)(?D a h**)(t) + ^(V„<" .(/i ) (n<)i 
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where we have omitted unimportant numerical factors and used the background Minkowski metric (|75|) to write the 
result. As usual, the notation f(t) stands for f(t, r = 0) for any function on spacetime. In ([55]) the quantity D a 
stands for the components of the derivative operator that acts on spacetime indices as the metric compatible one and 
on the internal indices as the covariant derivative -0,4b ■ To write the last term in ([85]) we have used the relation (|86p . 
Further, the covariant derivative D was replaced in this term by the usual metric compatible one because there are 
no internal indices to act on. A detailed derivation of (|89|) is given in the Appendix. Now using the property (188|) we 
see that the derivative operator D a in D a h l; > can be replaced by the partial derivative one (in place of which we 
can also use the metric compatible derivative operator as there are no spacetime indices in this quantity to act on). 
So, we get our final result for the self-dual term: 

J dt Qg^xrv^'xt) + ^(v„n(i)) • (so) 

Note that in the quantity Vah 1 ^ one first evaluates the derivative and only then evaluates the result at r = 0. 

Let us now treat the second, T* sd term in (|83[) . As we have already mentioned, this term is exactly the same as it 
is in GR. This is demonstrated in detail in the Appendix. So, we have for it: 

J dt (m(t) u a u b - im(t) g ab j (V a £ b ) (*). (91) 

Here u a is the vector tangent to 7, and the quantity in the brackets is simply the traceless part of the standard particle 
stress-energy tensor T ab = mu a u b . 
Now adding ([30]) and (JHU) we get: 

J dt (hjv(t){eV a W){t) + (m(t)u a )(u b V b ^)(t^j = 0, (92) 

where we have rewritten the second term in a suggestive form. Integrating in this second term by parts, and using 
the fact that £ a is an arbitrary vector field we get: 

u b V b {m(t)u a ) = ^Q ij (t){V a h ij ){t). (93) 

Equation (|93|) is our main evolution equation for "small bodies". It is instructive to see how the GR case gets 
reproduced. In that case there exists a unique background metric such that = <5 y . The right-hand-side of |93|) 
then vanishes and we get the usual u b W b (m(t)u a ) = 0, which implies that m(t) is constant along the worldline of the 
"small" body, and that this worldline is a geodesic. 

As a check of our result (|93|) we must make sure that it is invariant under conformal transformations of the metric 
used two write it. Indeed, in the theory under study only the conformal class of metrics is well-defined, not the 
metric itself. To see that the equation (|9"3")l is conformally invariant we recall the transformation properties of all 
the quantities: m — > fi — 1 m, Q lJ — > f2 3 Q y ', /i y — > fl^^h^ and finally u a — > ilu a , the last one following from the 
normalization condition g ah u a u b = 1. Thus, the quantity mu a is conformally invariant, and we only need to worry 
about the change of the metric-compatible derivative operator. For brevity we drop the argument indicating time 
dependence in all the formulae and get: 

Vf,(mit a ) — » Vb(mu a ) — TO(uf,V a lnfi + u a V b In £1 — g a bU c V c lnil). (94) 

This means that the left hand side of (|9"3"|) transforms as: 

u h V b {mu a ) -> rrVv^m-Ua) - fi _1 mV In ft. (95) 

On the other hand, the right-hand-side of (|93p transforms as: 

- lQ ij v a h ij -> -in 3 Q«v„jr 4 /^ = (96) 

which shows that the equation is indeed conformally invariant. 
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VII. INTERPRETATION 



In the previous Section we have obtained the "small body" evolution equation for our gravity theory. The only 
difference with the GR case stems from the fact that a "small body" is allowed to have a non-trivial "internal stress" 
tensor Q %s '(t) which then interacts non-trivially with the "non-metric" part of the background. 

From (j93|) we see that, because of the non-zero right-hand-side, the motion does not seem to be geodesic. Interest- 
ingly, when the background is metric, i.e. there exists a metric in which h lJ = 8 lJ , the evolution is geodesic even in 
case the body has non-trivial "internal stress" Q lJ (t) . Thus, it is only when the two-form field background in which 
the body moves is "non-metric" that we get (apparent, see below) deviations from geodesic motion. 

To give an interpretation to (|93|) . let us multiply this equation by u a . Due to the normalization condition satisfied 
by this vector field, we have u a \7bu a = 0, and so: 

u a W a m(t) = ^(t)(u a W a h^)(t). (97) 

Thus, the conservation equation no longer implies that the "mass" of the body is constant along the trajectory, but 
tells us something different. To see what, we will need to make an additional assumption about the nature of the 
"internal" stress-momentum of our body. Thus, we shall assume that the tensor Q^(t) that appears in (|97[) is 

independent of the direction of motion of the body. That is, we assume that the self-dual part of the body's stress- 
momentum two-form is (l/2)Q l; > with the same Q u no matter along which trajectory the body travels. Note that 
in general relativity this is true, with the self-dual part of T l being given by T' l sd = (m/6)B l . Our assumption may be 
motivated by considering how the self-dual part of T % arises from some matter action via (|5U|) . Indeed, T l sd arises from 
a term of the form (1 / '2)Q y B l A B^ , where Q is some matrix possibly depending on the "non-metric" components 
of B % . This suggests that the self-dual, or "internal" part of the stress-energy two- form should only depend on the 
internal composition of the particle, and not on its motion. 

If one makes this well- motivated assumption, then (|9T[) must hold for any choice of the vector u a . Thus, the 
following equation must hold: 

dm = ^Q ij dh ij . (98) 

Let us now recall that we have seen a similar equation before, equation (|32[) . in Section HTT1 that dealt with the vacuum 
modified gravity. There it implied that the tensor $ u that arises in the decomposition of the curvature into its self- 
and anti-sclf-dual parts must be given by a derivative of some potential function with respect to the internal metric. 
We see that a similar relation must be true here: 

«*« = 4 (tsT) «>■ 

Thus, the evolution equation implies that for each body there must exist a function m{h i: >) of the internal metric, 
such that the tensor Q lJ describing the self-dual part of the stress-momentum of this body is given by the partial 
derivatives of the mass function with respect to the components of the internal metric. Because of this, we shall no 
longer write the argument indicating the time dependence next to to, as we now interpret the mass as a function of 
h 1 ^ , which later must be evaluated on (t) . 

Let us now consider the "mass" function m(h l i) to be a function of the trace Tr(h) of the internal metric and the 
traceless part H % i . Then, using (f3"S"]) . as well as the fact that m(h) = Tr(Qh) we see that the function m(Tr(/i), H) is 
a homogeneous function of degree 1/4 in Tr(h). Therefore, we can write: 

(Tr(/i)W4 

m(h) = fh y ' J — W(H), (100) 

where W(H) = W(H t; i) is a dimensionless function of the traceless part H % i of the "internal" metric h 1 ^ normalized 
as W(Q) = 1, and to is a quantity with the dimensions of mass. Note that the formula (|100p is consistent with 
the transformation properties to — > f2 _1 m, h — * f2~ 4 /i of the quantities under the conformal transformations of the 
background metric. It should also be compared with an analogous formula (|38[) on the gravity side. 

Let us now see what the fact implies for the motion of the body. We can now replace the right-hand-side in 
(|93| by V a rn to get: 

u b V b m(t)u a = (V„m)(t). (101) 
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Recall that this equation is conformally invariant, with the mass function transforming as m — > f2 _1 m. Note also 
that the mass function is now defined not only along the trajectory, but everywhere, and we can use the conformal 
freedom in choosing the background metric to select a special metric in which m = in is a constant. Then in this 
metric, whose conformal factor is defined by the condition 

TT{h){W{H)f = 3 (102) 

the body moves along spacetime geodesies. Note that when the background is metric H = the mass rh is the usual 
mass of the particle as we know it in general relativity. So, similar to what we saw happening in the case of pure 
gravity, the departure from the familiar behaviour is parametrized by a single dimensionless function W(H) of the 
traceless part of the internal metric. 

It remains to discuss an interpretation of the function W(H). This function is set by the coupling of the matter 
component in question to the two-form field. In principle, it can be arbitrary, and moreover different for different 
matter components. In previous studies, see e.g. [10|, of the theory it was shown that the non-metricity H lJ of the 
two-form field background is controlled by the Weyl curvature, and that we can expect this non-metricity in the Solar 
system with its weak curvatures to be extremely small. This means that even if the function W{H) was different 
for different species of elementary particles one would not notice this in the Solar system. However, having different 
W(H) for different types of particles would mean that the weak equivalence principle was violated in the theory. 
Indeed, in this case each type of particle would travel along geodesies of its own metric and the universality of free 
fall would not hold (in regions of non-metricity). While this is an interesting theoretical possibility naturally provided 
by our theory, it is much safer to require the theory to respect the weak equivalence principle and postulate that the 
function W(H) is universal and same for all particles and composite bodies. In our further development of the theory 
we shall assume this to be the case. 

A universal function W(H) leads to one further interesting consideration. Indeed, the assumption of the universality 
of W(H) is the assumption that the dependence of all the mass terms in the matter Lagrangian on the H v components 
of the two- form field is the same. Let us now imagine that we have coupled matter not to gravity as in (|63|) , but to 
the topological BF theory: 

S top [B, A, 4> m ] = J B l AF l ~ 4nGS m [B, 4> m ]. (103) 

Now assuming all our matter fields <f> m to be quantum, one may take the vacuum expectation value of the last material 
term in the Lagrangian to get a Lagrangian that only depends on the two- form and the connection fields: 

S eff [B,A] = J B l AF-4iG(5 m [B,W). (104) 

The "quantum average" of the matter action must be an integral of a four-form that can only be built from the 4-forms 
B l A B 3 . However, these 4-forms are proportional to the internal metric h tJ and the volume form B % A B % . Thus, 
the four-form in question must be proportional to the volume form. Further, since all the mass terms in S f m [_B,0 m ] 
depend non-trivially on the components H lJ of B, the quantum average will depend on the function W(H). Thus, 
we see that the quantum average in (|104p must be equal to: 

S eff [B,A] = J B i AF' -^llu eff {H)(B l AS 4 ), (105) 

where A e ff is the effective cosmological constant, and U e ff(H) is some effective potential normalized so that £/ e //(0) = 
1. Here both A e ff and U e ff(H) depend on details of the matter Lagrangian and the form of the coupling W{H). 
What we have reproduced via this heuristic argument is precisely the gravitational action (|41|) with the potential 
U e f f (H) and the cosmological constant A e f / . This strongly suggests that the potential function U (H) appearing in 
gravity must be related to the mass function W(H) appearing in the matter sector. However, to find this relation one 
must perform a complicated quantum computation. In this paper we shall treat both functions as phenomenological, 
but one should keep in mind that it should be possible to relate them in the final theory. 

VIII. COUPLING TO THE GENERAL STRESS-ENERGY TENSOR 

Considerations of previous Sections fixed the form of the stress-momentum two-form of a "small body" . However, 
in order to be able to develop physical consequences of our gravity theory it is necessary to describe how arbitrary 
types of matter couple to it. Fortunately, the above considerations on the form of the stress- momentum two-form of 
a "small body" allow us to describe coupling to generic matter. 
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Since our theory respects the weak equivalence principle and there is a preferred metric in which test bodies 
move along geodesies and any metric is locally flat, it is natural to postulate, as in general relativity, that all non- 
gravitational physics in our theory is the same as in flat space. In particular, it is natural to assume that the effect of 
any matter on gravity is still characterized just by the usual stress-energy tensor of matter. Of course, the coupling 
of this stress-energy tensor to gravity may and will be different in the theory under study. 

The question thus reduces to that of describing how the usual stress-energy tensor of matter couples to our gravity 
theory. To answer it, the following formalism will be useful. As we have already done in Section [Vl given a general 
two- form field background B l , it will be convenient to introduce a set of "metric" two- forms. To this end, let us 
choose a representative in the conformal class defined by and choose a tetrad 9 1 ,1 — 0,1,2,3. From the tetrad 
one can construct the two forms B IJ := 9 1 A 9 and take the self-dual part with respect to the indices IJ . Let us refer 
to the so(3)-valued two-forms obtains this way as metric. They are the two-forms of Plebanski formulation of general 
relativity reviewed in Section [Til As before we denote these metric forms by a bold letter. Thus, we get a two-form 
field B 1 satisfying the metricity condition: B^-ABi ~ 5^-, as well as the reality conditions, see ([7]). As before, we shall 
continue to use the underlined indices to refer to quantities taking values in the "metric" SO(3) bundle that we have 
introduced via B 1 . Now, given the metric forms B 1 defining the same notion of self-duality on two- forms as B l , the 
original two-form field B l can be represented as a linear combination of the metric ones: 

B i = A i 1B i (106) 

The quantities A" are defined up to SO(3) rotations and rescalings of the metric two-forms B 1 . Thus, the invariant 
information contained in them is that in 9 — 4 = 5 components, and we can parametrize a general two- form background 
B % by its metric two-forms B- and by the quantities A", modulo conformal and SO(3) transformations. We shall see 
that this parametrization is very convenient for practical computations. In particular, the internal metric h lJ is given 
in terms of the matrices A" by: 

h ij = A^A*. (107) 

Let us now consider the stress-energy two-form. The usual stress-energy tensor T a b can be decomposed into its 
trace Tr(T) and traceless T a i, — (l/4)g a ;,Tr(T) parts. As we have explained above, we would like the stress-energy 
two-form T % of our theory to be constructed from the quantities Tr(T) and T f, — (l/4)g a fcTr(T). It can be expected 
that the traceless part will enter into the anti-self-dual part of the stress-momentum two-form to be constructed, and 
the trace part will enter into the self-dual part. Our experience with the stress-momentum two- form of a small body 
suggests that the anti-self-dual part of T l is essentially unchanged, and is given by: 

Tl sd = (A- 1 )^, (108) 

where A" is the matrix introduced in (|106p . and T l , is the "metric" anti-self-dual stress- momentum two-form, see 
Section [IT] Indeed, we have seen in (jSTj) that for a small body, the anti-self-dual tensor A"T^ sd has the usual form 
of one in the metric theory. We assume that this generalizes to arbitrary matter, and later check that this choice 
is consistent with energy conservation. This solves the problem of coupling the traceless part of T a b to our gravity 
theory. 

The coupling of the trace part Tr(T) is also suggested by what happens in the "small body" case. Thus, we saw 
that for a "small body" 

/-',/ ^"/>"- (109) 

and that the internal stress-momentum tensor Q iJ is given by a derivative of a "mass function" (|100[) with respect to 
the internal metric. We shall keep the same relation in the general case and write: 

Q' J =Tr(T)^, (110) 

where Tr(T) is the trace of the stress-energy tensor (with dimensions of energy density), and the dimensionless function 
R m (h) is given by 

R m (h) = ^-U m (H), (111) 

where U m (H) := (W(H)) 4 is the matter sector potential. Unlike in (|100p . which uses (Tr(/i)) 1 / 4 W / (/f), we have 
now used the fourth power of this combination. This is necessary for the internal stress-momentum tensor Q 13 to be 
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invariant under conformal transformations of the background metric. With this choice of T(h) we get: 

Q ij = ^(T)(j^ + ~5VU m y (112) 

where, as before, U m is the Legendre transform (jl3")) of the matter potential U m (H). Using the identity similar to 
(|3"T|) . it is now easy to see that 

h ij dQ ij = Rm(h)dTr(T). (113) 

We will need this identity below when we discuss the energy conservation. 

The expressions (jl08[) . (|109l) and (|112p determine the coupling of a general stress-energy tensor T a t, to our gravity 
theory. The only extra input that needs to be specified on top of what is already present in general relativity is two 
dimcnsionless potentials U(H), U m (H) of the traceless part of the internal metric. The potential U(H) determines 
the dynamics of the vacuum gravity, and the material potential U m {H) is necessary to specify the coupling to matter. 

It remains to check that the coupling specified is consistent with the standard energy conservation. The only thing 
that needs to be verified is that there are no changes in the self-dual part of the stress-energy two-form. The anti-self 
dual part does not change. A detailed argument involves spinors and is given in the Appendix. In our theory the 
conservation equation for T l is given by Using (|109[) . the self-dual part of this equation becomes: 

^B l A D AB Q»Bi = ^(B l A &)D A ,Q ij , (114) 

where we have again used the identity (f2"9"|) . Passing to the description of the two- form field B l in terms of the metric 
two-forms B-, and taking note of the definition (|107[) of the internal metric we can write this as: 

it c (B*AB*)i/i«dQ« = ^ € (B^AB^)i? m (/i)d(Tr(T)/4). (115) 

where we have used the second equation in (f64|) to replace the covariant derivative Da b by the usual one, and used 
(|113| to arrive at the final expression. The only modification here as compared to GR is the presence of the factor 
R m (h) in this expression. We see, therefore, that the conservation equation holds in the metric in which 

Tr(h)U m (H) =3, (116) 

in agreement with our finding (|102|) in the previous Section. This establishes that the standard stress-energy tensor 
T a b conservation equation is consistent with the conservation equation (|64p for T l when the stress-momentum two- 
form is constructed from the components of T a \, as specified in equations (|108| . (|109|1 and ()112|1 . A more thorough 
discussion of the energy conservation (including a treatment of the anti-self-dual part) may be found in the Appendix. 

This finishes the formal development of our theory. We now have all the ingredients to study physics with it, as we 
know how to describe pure gravity, know how it is influenced by a general stress-energy tensor, and know how test 
matter moves in a given gravitational background. At the end, our gravity theory has all the standard ingredients 
of general relativity. The only new quantities that we have introduced are two dimensionless "potential" functions 
U(H), U m (H) depending on the components H 1 ^ of the two- form field. Choosing the metric conformal factor so that 
the condition (|116[) is satisfied gives us completely standard physics for matter fields. The only thing that changes is 
the coupling of matter to gravity, as well as the gravitational dynamics. 



IX. RECIPE 



We finish our exposition of the new theory with an "explicitly metric" formulation that is useful for practical 
computations. The reader, however, will not see a metric below, only two-forms B- constructed from the metric, 
similar to what happens in Plebanski formulation of general relativity reviewed in Section [TTJ A more conventional 
formulation of the theory that uses the spacetime metric explicitly is also possible, see [ll[ for a recent treatment. 

As in GR in Plebanski formulation, see Section [Til one starts with a metric and the corresponding set of self-dual 
metric two-forms, which we denote by B 1 . One then forms a general linear combination of the metric forms: 

B i = A « B i ( 117 ) 

introduces the "internal" metric 
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and finds its traceless part: 



Tr(h) 

The theory is specified by two dimensionless potential functions normalized as 



W- 7 

H^ = ——-S^. (119) 



U(H) = l + ^Tr(H 2 )+0(H 3 ), U m {H) = 1 + ^Tr(H 2 ) + 0(H 3 ), (120) 

where a, f3 are dimensionless parameters. It can be seen that, apart from those already available in GR, these are the 
only parameters that are of relevance for the linearized theory. The matrices A" are required to satisfy: 

Tr{h)U m {H) = 3. (121) 

After this is done, one finds the two-form field compatible connection Ab such that Da b B 1 — (note that the 
derivative operator Da only acts on the non- underlined indices). The field equations then read: 

A^F l (A AB ) = A ( A'--^— A-'i + -h ij U) BL + 2nGT (.V--^-^ A'i + -h ij U m ) Bl - 2irGTHBi, (122) 



3 J V 3 

where U,U m are the Legendre transforms (|43|) of the potentials U(H),U m (H), A is the cosmological constant, T is 
the trace of the standard metric stress-energy tensor of matter, T 1 ^ are the quantities constructed from the traceless 
part of the standard stress-energy tensor, and B 1 are the anti-self-dual metric two-forms. For example, for the ideal 
fluid: 

T = (p-3P), T« = (p + P)S ij +ie^ k u k , (123) 

where p, P are the energy and pressure densities and u l is (related to) the momentum vector. The limit to general 
relativity is obtained by making the gravitational potential infinitely steep, i.e. by sending a — > oo,H i:i — > so that 
the product aH %3 remains finite. When H 13 = the matrix A" is an arbitrary SO(3) one, for example the identity 
matrix, and it is evident that (|122p reproduces Plebanski equations. The only new ingredient in (|122p . apart from 
the usual metric and stress-energy tensors, are the quantities A" that change the gravitational dynamics and the 
coupling to matter. They are, however, non-dynamical, and their only job is to "twist" the theory as prescribed by 
two potentials U(H), U m (H). 

The discussion of the previous sections has demonstrated that the stress-energy tensor is still conserved in this 
theory in the usual way V a T a f, = 0, where V a is the metric-compatible derivative, and that test bodies move along 
geodesies. In both cases the relevant metric is the one that is used in the construction of the metric two-forms B-. 
One can use the formulae given in this section as a definition of the theory. A reader who finds this definition a bit 
contrived should consult earlier sections for a simpler, but more abstract description. The physical exploration of this 
theory is left to future publications. 

For applications it is sometimes more convenient to work not with the internal metric h 13 , but with the matrices A l - 
introduced above in (|117|) . Thus, let us describe an equivalent formulation of the theory in which the internal metric 
never appears and one works directly with the quantities A". Again, we start with a metric and the corresponding 
set of self-dual metric two-forms, which satisfy B 1 A B^ ~ S 1 ^. As before, we form a general linear combination 
of the metric forms (|117p . However, now instead of introducing two potentials U(H),U m (H), let us work directly 
with the combinations: R(h) = (Tr(h)/3)U(H), R m (h) = (Tr(/i) /3)U m (H), which are two (arbitrary) SO(3)-invariant 
functions of the internal metric h lJ normalized so that R(6) — R m (S) — 1 and homogeneous of degree one in the 
quantity Tt(h). Let us view these functions as those of A 1 -. Then they are arbitrary (normalized) functions of 
matrices A" that are invariant under left and right action of SO(3), and transform as R — > fl~ A R,R m — > Vl~ A R m 
when A - — > £1~ 2 A". In order for the metric used to construct B- to be the physical one, in which matter moves along 
geodesies, the quantities A" are required to satisfy the conditions: 

R m (A) = 1. (124) 

As before, one now finds the incompatible connection Ab such that Da b B = (note that the derivative operator 
Da only acts on the non-underlined indices). After this is done, one computes the curvature of the connection Ab- 
From (|122| we see that the right-hand side of the field equations contains matrices of the type df(h)/dh lJ . Let us 
note an identity: 
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where / on the right-hand-side is considered to be a function of the matrix A". Using this identity we can rewrite 
the field equations (|122|1 in a way that uses directly the quantities A". One gets: 

A^-FUAab) = a( -^-A'l) Bl+ 2ttGT ( -^^A l i] W- - 27rGT«Bi. (126) 

The interpretation of the quantities A, T and is as before. Thus, for the ideal fluid we have (|123|) . The formulation 
that works directly with the internal metric "triads" A" thus leads to more compact field equations as compared to 
(I122|) and may be preferable for some purposes. The "twisting" role of the scalars A" is particularly clear in the 
formulation (|126p . 

From the described Plebanski-like formulation (|126|) it may seem that the obtained field equations have little to 
do with the objects one usually works with in gravity, namely the spacetime metric and the stress-energy tensor of 
some matter that couples to this metric. However, we would like to stress that in the final formulation our theory is 
completely standard and works exactly with the same quantities. Thus, we have the physical metric and the matter 
couples to it in a completely standard way. The matter moves along geodesies of this physical metric and has the 
usual stress-energy tensor. What is non-standard is how field equations for this metric are obtained. To this end one 
introduces certain extra scalar fields A" and deforms Einstein equations in a way that does not generate any kinetic 
term for the scalar fields and is consistent with energy conservation. 

Our theory in its final form may be compared to the "modified source gravity" of [l2j , where the authors, following 
[l3T |. introduce a scalar field ip and consider a gravity theory described by the following simple Lagrangian: 

J^{R-U{iP))+S m [eWg,ct> m ]. (127) 

M 

Here g is a dynamical metric, but note that the matter couples not to g but to a conformally related metric e 2 ^g 
instead. Importantly, the field ip is non-dynamical, i.e. does not have a kinetic term. In vacuum, the theory is just 
GR with a cosmological constant. But in general the theory for the physical metric e 2 ^ is different from GR and, in 
particular, the stress-energy tensor of matter sources the Einstein tensor of e 2 ^g in a modified way. 

The theory we have considered in this paper is similar to (|127[) in the sense that the coupling of the stress-energy 
tensor of matter to gravity is modified. The modification also arises from introducing non-dynamical scalar fields, 
even though there is now a multiplet A" of them instead of a single one in (|127|) . However, unlike in the case of (|127|) . 
the pure gravity theory is modified as well, with this modification being controlled by the potential R(h). Another 
important difference is that, unlike in the theory (|127[) that modifies the homogeneous isotropic Universe solution, in 
our theory the scalars A" are set to SO(3) matrices in this case by symmetry, so the homogeneous isotropic cosmology 
is unmodified. But the principle according to which ()127[) is constructed is quite similar to that used in our theory. 
This is made especially clear by a recent reformulation [ll| of the theory that works directly with the spacetime 
metric. 

Let us finish this section with two more remarks. As we see from the field equations (|126[) . in general, the self- 
and anti-self-dual parts of the stress-momentum two- form, or, in other words, the trace and the tracefree parts of the 
stress-energy tensor of matter tensor appear on the right-hand side of field equations on a different footing. Indeed, 
there is an extra matrix multiplying the self-dual part proportional to T in (| 1 26[) . An interesting question is if there is 
any metric among the conformal class defined by the original two-form field B l such that both parts of its stress-energy 
tensor appear in (|126p in the same way. Thus, we are looking for a function R m {A) such that dR m /dA ~ A -1 , such 
that Rm(d) = 1 and which transforms under A — > f2~ 2 A as R m — > £!~ 4 i? m . This function is: 

R m (A) = (det(A)) 2 / 3 , (128) 

or, in terms of the internal metric R m {h) = (det(ft,)) 1 / 3 . With this choice of the matter side potential the field 
equations take the form: 

A"-F l (A AB ) = A Bl + AkGT\ (129) 

with T 1 given by its usual expression in Plebanski theory, see e.g. (fT"9|) for the case of the ideal fluid. Let us finally 
note that the condition (|124[) in this case is simply det(h) — 1, which defines the so-called Urbantke metric flij . Thus, 
we can rephrase the above discussion by saying that the Urbantke metric is distinguished in our theory by the fact 
that when matter couples to this metric the field equations take a particularly simple form (|129|) . Urbantke metric 
has recently played a distinguished role in a reformulation of this theory proposed in [ll| . However, unlike in this 
reference, instead of fixing the metric to which matter couples from the outset, we prefer to allow matter to couple to 
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an arbitrary metric in the conformal class of B l , and control this coupling by the matter side potential R m (A). It can 
be seen that a non-trivial coupling with R m (A) different from (|128p leads to new interesting physical effects absent 
in the case (|128j) . which in our opinion serves as a sufficient motivation to allow such a more general coupling. 

It is also interesting to note that one can obtain a simple but still non-trivial theory with both potentials fixed by 
taking, in addition to (|128[) . the gravitational potential to be given by the same expression R(h) = (dot (ft,)) 1 / 3 . The 
obtained theory has no adjustable parameters and its field equations read: 

A i -F i (A\b) = -AB- + 4nGT-. (130) 
3 

If not for the presence of the quantities A" the vacuum T 1 = version of these equations would be just the constant 
curvature condition. The presence of the extra scalars twisting this equation makes them more interesting, and, in 
particular, makes a non-trivial spherically-symmetric solution possible. However, in view of the fact that the modified 
gravity theory (|130|) does not have adjustable parameters it is likely to be in gross conflict with the standard gravity 
tests. Thus, the particularly simple version of the theory with both potentials fixes to be R(h) = R m {h) = (det(fe)) 1 / 3 
is likely to be of only academic interest. 

X. DISCUSSION 

With most of the discussion being embedded in the main text we shall only make some "philosophical" remarks 
here. The role of the metric in general relativity is two-fold. First, a metric defines the spacetime causal structure 
(lightcones at every point). However, to define the causal structure one only needs the conformal metric, i.e. the 
metric modulo conformal rescalings. Second, every spacetime point in general relativity is effectively equipped with 
a set of rulers and a clock. It is for this purpose of measuring spacetime intervals that one needs a metric per se, not 
just a conformal metric. While the propagation of light is a very basic process that can arguably make sense to be 
built into the very definition of the spacetime structure, the availability of rulers and clocks at every point is on a very 
different footing. Indeed, a measurement of distances and time intervals is a complex physical process that requires 
in each case a macroscopic physical system - a solid body or a clock. The very availability of rulers and clocks at 
every spacetime point is quite striking, for even the empty space, which is by definition void of anything material, is 
endowed in GR with this structure. In our opinion this is a very anti-Machian feature of GR: according to Mach's 
ideas only a relative description of material bodies in the Universe is possible, and an "empty" Universe filled with 
clocks should be approached with suspicion. One can argue that the notion of spacetime distance is a macroscopic 
one, and has no place in any reasonable microscopic description. The fact that macroscopic bodies behave as if to 
register spacetime intervals needs to be explained, not postulated. This is not so for the causal structure, as the tiny 
quanta of electromagnetic field constantly popping out of the vacuum can be argued to define the causal structure of 
spacetime by their very existence. 

It thus seems reasonable to try to formulate a theory of the gravitational field which is based on the spacetime 
conformal structure, not on the spacetime metric. Experimentally we only know that the gravitational field is a 
universal long-range interaction. As such it should be possible to think about it as occurring due to exchange of 
some massless particles which can thus have only two possible polarizations. In GR these two polarizations arise from 
the gravitational field of a spacetime metric with its ten components and an additional 4-parameter group of gauge 
symmetries - diffeomorphisms. However, it seems impossible to build a similar description on just the conformal 
structure, as it is specified by 9 components, which is an odd number, and so no straightforward scheme with gauge 
symmetries (which reduces the number of DOF by an even number) can bring 9 total components down to two 
physical. A theory of gravity that is based on just the conformal structure would also have the problem that there 
would be no preferred scales in it, so the world described by it would not be realistic. 

The theory that we have formulated in this paper describes spacetime geometry by specifying its conformal structure. 
The way this happens is that in addition to the conformal structure there are other fields in the theory - other 
components of the gravitational field. The total number of "components" of our basic two-form field is 18 - an even 
number and taking into account all the arising constraints it can be seen that the number of the arising polarizations 
of the graviton is still two. The vacuum theory is specified by one arbitrary dimensionless scalar function U(H) of 
a traceless symmetric 3x3 matrix H % K When one couples the theory to matter one has to introduce yet another 
arbitrary function U m (H) with similar properties. It is this matter sector potential function U m {H) that can be 
shown to supply the conformal factor that defines the spacetime metric in which test bodies move along geodesies. 
The limit to general relativity is obtained by setting H %3 = 0. 

As we have discussed, the gravity and material sector functions should be related, since at least in principle it 
should be possible to compute the gravity potential U(H) as induced by quantum effects involving matter. Indeed, 
in Section IVIII we have seen how a version of "induced gravity" scenario is possible in our setting. With the current 
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state of the development of the theory such a quantum computation remains beyond our abilities and there is no 
choice but to treat the two functions as two independent phenomenological parameters (or rather two infinite set 
of parameters) of the theory. It should be noted, however, that the linearized theory is only sensitive to the two 
leading parameters from the arising infinite towers of coefficients. A work analyzing the effect of modification on the 
cosmological perturbation theory is currently in preparation. It should also be noted that the spherically-symmetric 
solution of the described gravity theory is known, see e.g. [l(|. Implications of the modification for the motion of test 
bodies in the spherically-symmetric background need to be analyzed in light of findings of this paper. 

Let us conclude by expressing our amazement at how far it was possible to develop our modified gravity theory just 
by following its internal logic. Indeed, the theory started its life in Q as a rather complicated modification of the pure 
connection formulation of GR. However, as we have seen from the constructions presented in this paper, the theory 
turned out to be a very natural generalization of Plebanski gravity tightly constrained by "Bianchi" identities, both 
the vacuum theory and its coupling to matter. We have also seen that after specifying the matter sector potential 
U m (H) (or, equivalently, R m {h)) and thus specyfying to which metric in the conformal class defined by B % the matter 
couples, the theory takes an entirely standard metric form. 

It should also be emphasized how striking the results described in this paper are from the more familiar perspective 
of metric-based gravity theories. Indeed, it is commonly believed that in order to modify gravity one needs to 
introduce new propagating degrees of freedom. In addition, it is often said that the gravitational coupling to matter is 
constrained by the energy conservation, which leaves essentially no freedom. The modification described in this paper 
changes GR without adding to it any extra propagating modes. In addition, a similar modification of the coupling of 
the stress-energy tensor to gravity becomes possible, without any contradiction to energy conservation. 

In spite of these exciting results much more remains to be done. In particular, the question of coupling of fermionic 
matter directly to the two-form field remains open. As we have seen in this paper, it is not necessary to answer 
this question to obtain physical predictions of the theory, but it will certainly be an essential question when the 
quantization is attempted. We thus hope that the theory that started its life almost 20 years ago in Q will continue 
to be a source of interesting results for some more years to come. 

Appendix: Spinor techniques and energy conservation 

Let us, as before, select an arbitrary metric in the conformal class of metrics determined by B l , construct a tetrad, 
and then use it to identify the space of rank 2 mixed primed-unprimed spinors Xaa 1 with spacetime vectors A° (and 
also, using the metric, with spacetime one- forms A Q ). Thus, all spacetime indices are converted into spinor ones. In 
these notations our basic two-form field B % is a self-dual two-form B AB eA'B' (there is no component proportional 
to €ab which would correspond to the anti-self-dual part). Thus, the two-form field B l becomes described in this 
language by the metric that it defines, as well as by the quantity B AB , symmetric in the unprimed spinor indices AB. 
In the GR case B\ B B l CD ~ £A(c e \B\D)i but in general these quantities are arbitrary. The stress-momentum two-form 
T % is described by its self-dual T AB and anti-self-dual T A , B , components. 

The conditions (|83|) whose consequences we need to explore in spinor notations become: 

J dt (-T iAC D A 't*B AB + T iA ' c ' D A ^B\ B ) = 0, (131) 

where ^ AA is the spinorial representation of the vector field £ a , and Daa' is that of the covariant derivative operator 
Da b - 

Let us now use the form (|8ip of the stress-momentum two form. The anti-self-dual component of T l is given by: 

Tk'B' = {B-^ab^a-b; (132) 

where Taba'B' has the same form as in GR TabA'B' = mu(A\A'\ u B)B' > and is just the traceless part T a b — (l/A)Tg a b 
of the stress-energy tensor T a b of ideal pressureless fluid T a t — mu a Ub,u a u a = 1. The quantities B\ B is what A lJ 
become in the spinor notations, and {B~ 1 ) l AB is the inverse matrix satisfying: 

{B~ l ) iAB B AB = (B-^&cj, = e A{ ce\B\D)- (133) 

Let us thus look at the second term in (| 13 1[) . The compatibility equation DB l = takes in spinor notations the 
following form: 



D%B\ B = 0. 



(134) 
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This means that the quantity B AB can be taken out of the operator of the covariant derivative: 



rr^iA'C' j-> A c B r>i rpi A' C' r>i 7-. A <■ B rniA'C r>i V7 A /■ B 

1 U C' { iA' B AB - 1 tS AB L) C' { iA' - 1 B AB \J c ,t, A ,, 



(135) 



where the last equality is due to the fact that there are no internal indices in the quantity that the covariant derivative 
operator acts on, and so the derivative operator can be replaced by the usual metric-compatible one. We now use the 
form (|132p of the anti-self-dual part of the stress-momentum two form to conclude that the second term in (|131[) is 
given by: 



T 



A'C'r-jA C B 



AB 



(136) 



where we have used the fact that Taba'B' — mu (A\A'\ u b)b* ■ 

We can now substitute (fB"l| into the first term self-dual term in the conservation equation (|131[) . and use the 
compatibility equation to take the quantity B 3 Ac under the operator of covariant derivative. The first term becomes: 



~QVD£'a'Bi AC B 



AB- 



(137) 



We now use the identity 



which in the spinor notations becomes: 



to rewrite (|137j) as: 



B% A B% = -U EF h^, 



-O^D A '^ = -l Q vtf,D£hV ^Dfa, 



(138) 



(139) 



where we have used the relation (|86|) . Now using the identity (|88|) we see that we can replace the covariant derivative 
in the first term here by the ordinary one. The covariant derivative in the second term can be replaced by the metric 
compatible one as the quantity it acts on does not have internal indices. 
Combining it all together we get for the equation (|131|) : 







(140) 



where we have rewritten the second term in a suggestive way. We now use: 



.a'.b' 



1 



A'B' 



l A U B 



to get: 



or, in the usual tensorial notations: 



Q = J dt(^Q l ^ b \7 b h^ +mu b {u a \7 a £_ b )^j = J dt Qq^V, 



bh 13 - u a V a (mu b ) £ 



(141) 



(142) 



(143) 



where to obtain the second equality we have integrated by parts in the second term. Since the vector field £ a is 
arbitrary we can conclude that: 



u a V a {mu b ) = U)'*V,,h' ! . 



(144) 



This finishes our proof of (f9"3"|) . 

Let us now consider a proof of energy conservation. For this we take the energy conservation equation in the form 
(|64l) . which in spinor notations becomes: 



i%B\ B {~D CA T i 



c 



D AC 'T^') 



0. 



(145) 
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Let us transform the anti-self-dual part first. We can use the compatibility equation (|134[) to take the quantity B AB 
under the operator of covariant derivative. We get for this term: 

d*D AB ' B\ B T B A ' = £%V ab 't a a b „ (146) 

which coincides with the usual expression in the metric theory. 

Let us now analyze the self-dual term. Here we can replace the self-dual stress-momentum two-form by its expression 
(|84p in terms of the tensor Q tJ , and take the two-form B l out of the operator of covariant derivative. We get for this 
first term: 

- l^B AB B> c A D CA 'QV, (147) 



where Q lJ is now given by: 



««-T^S», (14S, 



with T being the trace of the stress-energy tensor T a i,. Now using (|138[) we get: 

~ti,h»D/Q* = ^fl(ft)V/(T/4), (149) 

where we have used the identity (1113|) . Combining the two terms above, and writing the result (up to an overall minus 
sign) in usual vector notations we get: 

fl(fc)£°V„(T/4) + e Q V b (T ab - (l/4) 5ab T) - 0. (150) 

This coincides with the usual conservation equation V a T {, when the metric is chosen so that R(h) = 1. This finishes 
our demonstration of the fact that the usual energy conservation holds in our theory. 
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